Nonasymptotic Probability Bounds for Fading 
Channels Exploiting Dedekind Zeta Functions 

Camilla Hollanti, Member, IEEE, Emanuele Viterbo, Fellow, IEEE, and David Karpuk, Member, IEEE 



m 
o 

(N 

Oh 
< 
00 



CO 



(N 
> 

in 

m 



Abstract 

In this paper, new probability bounds are derived for algebraic lattice codes. This is done by using the Dedekind 
zeta functions of the algebraic number fields involved in the lattice constructions. In particular, it is shown how 
to upper bound the error performance of a finite constellation on a Rayleigh fading channel and the probability of 
an eavesdropper's correct decision in a wiretap channel. As a byproduct, an estimate of the number of elements 
with a certain algebraic norm within a finite hyper-cube is derived. While this type of estimates have been, to some 
extent, considered in algebraic number theory before, they are now brought into novel practice in the context of 
fading channel communications. Hence, the interest here is in small-dimensional lattices and finite constellations 
rather than in the asymptotic behavior. 

Index Terms 

Bounded height, Dedekind zeta function, lattices, norm forms, number fields, PEP, theta series, union bound, 
unit group, wiretap channel. 

I. Background and related work 
It has been well known for many years that number field lattice codes provide an efficient and robust means 
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^ . for many applications in wireless communications. We refer to d for a thorough introduction to the topic. More 
T"! 'recently, wiretap channels and number field based codes have been under study. Gaussian and fading wiretap 

> : 

channels have been considered for example in ||3l, 0], J5], JSJ, |7). In (8J the authors propose to use number field 
^ lattice codes, which will also form the basis for our study and constructions. This paper can be seen, on one hand, 
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as a continuation of (9], Q, where analysis on lattice codes in fast and block fading channels was carried out 
based on various explicit code constructions and, on the other hand, of iTTOl . ifTTI . where Vehkalahti and Lu showed 
how the unit group and diversity-multiplexing gain trade-off (DMT) of division algebra-based space-time codes are 
linked to each other through inverse determinant sums, and also demonstrated the connection to zeta functions and 
point counting. As continuation to ifTOl . IfTTI . the authors later showed that the density of unit group completely 
determines the growth of the inverse determinant sum, see lTT2l . |[T3l . 

Our work differs from lfl"3l in that it is targeted towards practical SNR region and small delay rather than to 
asymptotic behavior. In addition to the pairwise error probability case, we apply the same methods to the wiretap 
channel, where an eavesdropper is trying to intercept the data. While in lfT3l the authors concentrate on the number 
of units in a finite spherical subset of a lattice (as this is known to be the dominating factor), here we bound each 
individual term in the norm sum. This will enable us to estimate the number of points with each norm within a 
hyper-cubic constellation resulting in finer probability bounds. 

In summary, the contributions of this paper lie within 

• deriving (nonasymptotic) bounds for the probability expressions related to Rayleigh fading channels and to 
(Rayleigh fading) wiretap channels by using Dedekind zeta functions, 

• finding more accurate bounds through geometric analysis based on the unit lattice, Dedekind zeta functions, 
and bounded height norm sums, 

• deriving an estimate to the number of elements with certain norm as a byproduct, 

• demonstrating the accuracy of the estimate through numerical examples and showing that the estimation error 
is very small for small dimensions. 

While bounded height (cf. Def. [5]) estimates of the same type as the ones derived in this paper are known in 
algebraic number theory lTT4l . ITT31 . they are far from being standard or well-known. Hence, we hope that our 
derivation and then the practical use of such estimates will boost further research leading to yet tighter bounds. 

One should keep in mind that in general the estimation error is not negligible when the lattice dimension grows. 
However, as we apply these estimates in the situation where the lattice dimension corresponds to the decoding 
delay, the error is less severe. We will show that the estimate is very good when the dimension is relatively low 
and hence the delay short. Notice that the lattice dimension is not limiting the data rate as we can always increase 
the constellation size by choosing a bigger hyper-cube. Actually, the bigger the cube, the better our estimate will 
be, since the edge error effect becomes more negligible. Another limitation to the lattice dimension is forced upon 
by decoding - it is known that the complexity of any ML decoder such as a sphere decoder grows exponentially 
in the lattice dimension. We refer the reader to lTl4l . lTT~5l for more details regarding norm forms and the treatment 
of related error terms. We also point out the similarity of the constant term obtained here to the one in lfT~5l . which 
well demonstrates the accuracy of our derivation. 

The rest of the paper is organized as follows. In Section JTJ we provide some algebraic preliminaries related to 
number field lattice codes. Section JII] shortly introduces lattice coset coding employed in wiretap communications 
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and revisits the probability expression for the fast fading wiretap channel as well as for the typical Rayleigh fading 
channel, finally unifying the treatment of both expressions as one to simplify the computations in the rest of the 
paper. In Section |IVJ first bounds are derived using Dedekind zeta functions. We then refine the bounds in Section 
IVl where geometric analysis is carried out in order to estimate the number of constellation points with certain 
algebraic norm. The accuracy of the estimate in low dimensions is demonstrated in Section [VT] Conclusions are 
provided in Section I VII I 

II. Algebraic preliminaries 

Lattices will play a key role throughout the paper, so let us first recall the notion of a lattice. For our purposes, 
a lattice A is a discrete abelian subgroup of a real vector space, 

A = ZPi © Z/3 2 • • • © Zp t C R", 

where the elements /3i, . . . , (3 t G R n are linearly independent, i.e., form a lattice basis, and t < n is called the rank 
of the lattice. Here, we consider full (t = n) totally real lattices arising from algebraic number fields (see Def. [3] 
below). For more details on lattices and lattice codes, refer to 0. 
Definition 1: The minimum product distance of a lattice A is 

n 

rf Pl mm(A) = min 1 \\xi\, 

where x = (xi, . . . , x n ) G A. 

Definition 2: Let K be a number field. A real embedding is a field homomorphism a : K ^ R. A complex 
embedding is a field homomorphism a : K C such that a(K) C R. 

Definition 3: Let K/Q be a totally real number field extension of degree n and a±, . . . ,ak its embeddings to 
R. Let Ok denote the ring of integers in K. The canonical embedding tp : K R n defines a lattice A = iJj(Ok) 
in R n : 

ii(x) = (ai(x), . . . ,a k (x)) G ^(O k ) C R n , 

where x G Ok- 

Definition 4: The algebraic norm of x G K is defined as 

n 

i=l 

We abbreviate N(x) = N k /q(x) whenever there is no danger of confusion. 
If x G Ok, then Nk/q(x) G Z. Hence, we have that 

d p ,min{^{C>K)) = min \N k /q(x)\ = 1. 

In what follows, a cubic constellation will be used, bounding the size of the vector components in the canonical 
embedding. To this end, we define the height of an algebraic integer as follows. 
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Definition 5: The height of x G Ox is 

H(x) = max |<7j(a;)|. 

l<i<n 

We extend this notation to the height of a principal ideal X = (x), x E Or-, in a natural way by identifying an 
ideal X with its minimum-height generating element and simply defining 

H(X) = mm{H(y)\X = (y)}. 

Let us denote by (n,^) the signature of K, i.e., 

[K : Q] = n = ri + 2r 2 , 

where r\ is the number of real embeddings K ■— )■ R and r2 is the number of conjugate pairs of imaginary 
embeddings K C. The group 0^- of units of 0^ is described by the following well-known theorem, repeated 
here for the ease of reading. 

Theorem 1: (Q~6j Dirichlet Unit Theorem 1.9]) Let K be a number field and let (n,^) be the signature of K. 
There are units e±, . . . , e ri+r2 +i € such that 

0£ = Wk x (ei> x ••• x (e ri+r2 _i) 

where Wk is the group of roots of units in K. The ej are called a fundamental system of units for K. 

The fundamental units are used for defining the regulator of K. Let {ei, . . . , e r } be a fundamental system of 
units for K, where r = T\ + r% — 1. Consider a matrix 

4 = (log | a, (e^ | j ) 

for 1 < i < r and 1 < j < t\ + r2, and where we have used the notation 

{M if 1 < j < ri, 
|x| 2 if ri < j < ri + r<i. 

Here |x| is the usual absolute value on C, and o~\, . . . , a ri are all the real embeddings, while o> 1+ i, . . . , cr ri+r2 are 
a set of representatives of the imaginary embeddings. 

Definition 6: The regulator px is the absolute value of the determinant of any r x r minor of A. It is independent 
of the choice of the fundamental system of units. 

The volume of the fundamental parallelotope of the log-lattice Ai og generated by A is expressed in terms of the 
regulator as 

vo\(Ai og ) = p K Vn + r 2 . 
In the case of a totally real number field we have 

vol(A io9 ) = p K Vn- 



5 



The regulator is a positive real number that in essence tells us how dense the units are. The smaller the regulator, 
the denser the units. Regulators can be easily computed by the Sage computer software ifTTl . 

Let us conclude this section by defining the Dedekind zeta function and its truncated form. 

Definition 7: The Dedekind zeta function (cf. lfT6l p. 37]) of a field K is defined as 

where X runs through the nonzero integral ideals of Ok- The sum converges for > 1. Since N k /q(Ok) = 1, 
we always have 

Ck(s) > 1 

for all integer exponents s > 2. From now on, we assume 2 < s G Z since these are the interesting values for the 
applications under study in this paper. Namely, s = 3 corresponds to the wiretap case and s = 2 to the PEP case. 
Values s > 3 become relevant in the multiple-input multiple-output (MIMO) case. Convergence is guaranteed in 
all these cases. 

The Dedekind zeta function can be written as a Dirichlet series 

fe>i 

where a& = for those k that do not appear as a norm of an integral ideal. 

Definition 8: The bounded height Dedekind zeta function of a field K is defined as 

1 Nm a 

0<H{X)<m 7 fc=l 

where 

iV m = max{/c | k = N(l), H(X) < m}. 
Also the zeta functions as well as the Dirichlet coefficients can be computed by Sage IfTTl . 

III. Probability expressions and inverse norm sums 
A. Wiretap channel and the probability of Eve 's correct decision 

In a wiretap channel, Alice is transmitting confidential data to the intended receiver Bob over a fading channel, 
while an eavesdropper Eve tries to intercept the data received over another fading channel. The security is based on 
the assumption that Bob's SNR is sufficiently large compared to Eve's SNR. In addition, a coset coding strategy 
lfl"8l is employed in order to confuse Eve. 

In coset coding, random bits are transmitted in addition to the data bits. Let us denote the lattice intended to 
Bob by Aft, and by A e C A& the sublattice embedding the random bits. Now the transmitted codeword x is picked 
from a certain coset A e + c belonging to the disjoint union 

A b = uf =1 A e + cj 
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embedding k bits: 

x = r + ceA e + c, 

where r embeds the random bits, and c contains the data bits. 

We assume the fading is Rayleigh distributed and that both Bob and Eve have perfect channel state information 
(CSI-R), while Alice has none. The details of the channel model and related probability expressions can be found 
in 0. 

Next, let us recall the expression P c >e of the probability of a correct decision for Eve, when observing a lattice 
A e . For the fast fading case (U Sec. III-A], 

n/2 

0^xGA e 

where j e is the average SNR for Eve assumed sufficiently large so that Eve can perfectly decode A e . This is a 

reasonable assumption, as Eve is assumed to have perfect CSI. Here A& denotes the lattice intended to Bob, and 

A e C A;,. It can be concluded from ([3]) that the smaller the sum 

n 1 

e n uTj3' 

O^xGA, i=l 1 *' 

the more confusion Eve is experiencing. 

As a construction method, the authors of ||8] propose to use the canonical embedding of the ring of integers Ok 
(or a suitable proper ideal X C Ok) of a totally real number field K over Q. The field K is chosen totally real to 
achieve full diversity. More precisely, if x G Ok, the transmitted lattice vector in the fast fading case would be 

x = iP(x) = (a 1 (x),a 2 (x),...,a n (x)) G A e C R n , (4) 

where tp denotes the canonical embedding (cf. Def. [3]) and o-j are the (now all real) embeddings of K into R. 
The corresponding probability of Eve's correct decision © yields the following inverse norm power sw^ to be 
minimized O Sec. III-B]: 

Sk = ^2 Tat MiT' ( 5 ) 

Remark 1: The above sum Sk may not converge, since infinitely many elements can have the same norm. This 
happens e.g. when the unit group is infinite, which is the case for all field extensions other than the trivial one and 
imaginary quadratic fields. In practice, however, we always consider finite signaling alphabets, so the sum becomes 
truncated and converges regardless of the field. 

Thus, we conclude by the following natural definition of the truncated version of Sk, 

S «^= E |aw.)|3 = Eg> 

0<H(x)<m 1 A/l4V n fc=l 
'The behavior of the probability of Eve's correct decision in some example cases has been analyzed in (9). 
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where N m = maxjfc | k = N(x), H(x) < m} is now denned for element norms similarly as we did for ideal 
norms in the case of the truncated zeta function. 

Remember that H(x) was only defined for algebraic integers, so the above notation implicitly contains the fact 
that x e K - 

B. Rayleigh fading channel and the pairwise error probability 

Similar sums can be used to get bounds for the pairwise error probability for the fast fading channel when 
employing the same number field code as in (0]). Let us use the following Dirichlet series notation for the PEP (see 
e.g. El), 

P =A V - = ^V^ (7) 

e 7™ ^ \N K/0 (x)\ 2 7 "^fe 2 ' K) 

where 7 is the SNR and c is a fixed constant. We shortly denote the truncated version (ignoring constants) as 

-, N m , 

s K , PEP (m)= £ ,]vwiyF = E#- m 

0<H{x)<m 1 n k=l 

C. Unified treatment of the norm sums 
Let us use a more general notation 

Now we can carry out the analysis for both the sums © and © by simply denoting 

S K (s,m)= W M|. = Ep (10) 

0<H(x)<m ' K /^ X " k=l 

and then inserting s = 2 in the Rayleigh fading case, and s = 3 in the wiretap case. 

IV. First bounds 

In this section, we will derive lower and upper bounds for the inverse norm power sum (cf. (|6), ©) by using 
the Dedekind zeta functions (cf. Def. |7]). 

Similarly to the zeta function, we trivially have that Sk{s) > 1. Albeit straightforward, the following result gives 
us a nontrivial lower bound / 1 for the sum Sk(s). 

Proposition 1: (Lower Bound) 

Assume that Ok is a principal ideal domain (PID) and A e is as above with x G Ok- Let m denote the maximum 
height included in the sum. The Dedekind zeta function (k(s) (1 < s € Z) provides us with a lower bound for 
bounded the height sums S*x(s,m), i.e., 

S K (s,m) > ( K {s,m) > 1, 

If K is the field or rationals or an imaginary quadratic number field, also the unbounded sum will converge. 
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Proof: Note that N k /q(I) = [Ok : Z] G Z + , and that in the zeta function the summation only goes through 
the (integral) ideals of Ok, whereas in Sk we sum over all the algebraic integers of K. 
Let us denote by 

SK = Y. b fs and ^{s) = y^ys 

n>\ n>l 

the Dirichlet series |[T6l p. 31] of Sk and (k{s). Denote further by 

A = {k | a k / 0} C Z+ 
the set of values k that appear as norms in Ck( s ), an d by 

B = {k | b k ? 0} C Z+ 

the set of values that appear as norms in Sk- As K is a PID, we know that N k /q(Z) = mino^gi \N k /q(x)\, 
and that N k /q(Z) = N K /Q{(a)) = \N K /Q{a)\, where a is the generator of 1. Hence, we have that 

A = B. 

Further, we easily see that b k > ak- Namely, if k appears as a norm for distinct ideals Xj = (a,), i = 1, . . . , a k , 
it then appears as a norm at least for the (distinct) elements G Ok- In addition, we may have an element 
a G Ok, a. ^ ajVi with the same norm n. On the other hand, if k appears as a norm for distinct elements 
ai G Ok, i = 1) • • • , b/., it cannot appear as a norm for an ideal for more than a k = b k times, since Ok is a PID. 
This proves the claimed lower bound for the infinite sums. Now it is obvious that the bound also holds for the 
bounded height sums, where the maximum height is the same for both sums. 

■ 

Remark 2: According to Dirichlet's Unit Theorem (cf. Prop. Q3 the rank of the group of units in Ok is r = 
f \ + f2 — 1, where (ri,r2) is the signature of K, that is, n = r\ + 2r2- Therefore the infinite sum Sk does not 
converge for totally real number fields 2 Q ( r i > 1)> as tne y have an infinite unit group (cf. Remark [T). Hence, in 
the present application case where K is totally real, we indeed always need to consider truncated sums. In practice 
the constellations are always finite, hence this poses no real restrictions. 

Remark 3: Dedekind zeta functions have also been used in ifTOl for studying the DMT of the multiple-access 
channel (MAC). Interesting subsequent work has been carried out in ifTTIl . |[T3l . lTT2l . 

Next, let us denote by 

SK(m) = Y J h i 

k=l 

the truncated sum, where b' k = for those k that do not appear as a norm for x G Ok, H(x) < m. An upper 
bound for the truncated sum is achieved from the truncated Dedekind zeta function (^(s,m). 
Proposition 2: (Upper Bound) Let Ok be a PID. Then we have that 

SK(m) < max{6' fc | k < N m } ■ (k(s, m). 
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Proof: First, note that if = 0, then bk = and hence b' k = since Ok is a PID. Now a simple computation 
gives us 

S K (s,m)= Ys 

k<N m 

< m&x{b' k | k < N m } U 

b' k ^0,k<N m 

< max{b' k | k < N m } ■ ( K (s, m). 

In the second step, the summation is only taken over the values of k for which b' k ^ 0, i.e., k appears as a norm 
for some element x. Hence, for all the terms l/k s included in the sum / 0. In addition to this, the truncated 
zeta sum may contain other terms (for which 6^/0 but b' k = due to the energy limit) so we indeed get an upper 
bound. ■ 
Remark 4: We anticipate that the above bounds are not very tight, especially when extended to the non-PID 
case. Hence, our goal in the next section is to derive tighter bounds arising from geometric analysis. 

V. Probability bounds from geometric derivation 

Let us consider an algebraic lattice generated by the canonical embedding of the ring of integers Ok of a totally 
real algebraic number field K = Q(9) of degree n, 

A = {x = (xi,...,x n ) = 

= (<7i (a), (72(a),- • • ,a k (a))\a £ O k }, 

where the cr^'s are the n field homomorphisms. In the following we will indistinctly talk about lattice points and 
algebraic integers. 

We will consider a finite constellation S = A n B where the bounding region is a hypercube of side 2R centered 
at the origin. 

The algebraic norm of an algebraic integer is equal to the product distance of the corresponding lattice point 
from the origin, 

n n 

N(a) = Y[a i (a) = ~[[x i . 

i=l i=l 

Since the field norms of algebraic integers to Q are integers we observe that all the lattice points on the hyperbolic 
naps 

n 

| \Xi\ = k 
i=l 

correspond to algebraic integers with the same absolute norm. When k = 1 we have the group of units of K which 
are fully characterized by the Dirichelet theorem. 



10 



In order to evaluate the union bound on the error probability and the probability of the eavesdropper's correct 
decision, we now use a slightly different notation for the truncated sum as a product distance theta series, 

S,k>l 

where n k is the number of constellation points at product distance k from the origin when the sums are limited to 
the points in the finite constellation S. 

We need to count the number of lattice constellation points on each hyperbolic nap. We will consider the following 
logarithmic coordinate system defined by 

Xi = log | |, i = 1, . . . ,n. 

In this logarithmic space the hyperbolic naps become hyperplanes T-L k perpendicular to the all ones vector. In 
particular, the units are on the hyperplane passing through the origin of the log-space. We also note that the number 
of hyperbolic naps that are mapped to the same lattice is given by the number w of roots of unity that are in the 
field K. 

The bounding box B is transformed in the log-space to a semi-infinite rectangular region with a corner at 

(log(i?),...,log(i?)) 

denoted by log(£>) = (— oo, \og(R)] n . 
Since k is positive, the intersection 

S k = \og{B)C\U k 

is nonempty for all absolute norms 1 < k < R n . 

Dirichelet's Theorem proves that the units of K form a lattice Ai og in the hyperplane 

n 
i=l 

with a fundamental volume Vol(Ai og ) = pKy/ r i + r 2> where px is the regulator of the number field (cf. Def. [6) 
and n = r\ + 2r2 (n = r\ when K totally real). 

This suggests a simple and well-known way to estimate the number of units in the box as 

lt!VOl(5l) 
Tl\ = 

_vol(Ai og )J ' 

where |_-J denotes the floor operation. 

For non units (k > 1) the problem is more complicated. We need to count the number a k of principal ideals 
Ij{k) = (ctj) of norm k generated by the elements ay of absolute norm k. Since N(au) = N(a) for all units u 
and the norm of a principal ideal is equal to the absolute norm of its generator we can conclude that the set of 
points on the hyperplane S k ,k > 1 is a union of translates of the logarithmic lattice. Then we can estimate the 
number of points of norm k > 1 in the box as 

wa k vol(S k ) 
vol(Ai og ) 
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Definition 9: If we denote by bk the exact number of points of absolute norm k, then we define an error function 

fh(R) = LK-MJ, 

which is dependent on the edge length R. 

Remark 5: It is known that the error function grows quite large when the dimension of the lattice grows. We 
will illustrate the size of the error function later in Section ED We are not making any claims as to how tight our 
estimate n& is asymptotically, but rather wish to show that with small values of n with practical interest it will give 
us a good estimate. 

Let us now focus on the derivation of vol(<Sfc). We first note that Sk is the basis of a hyper pyramid Vk with 
vertex in (log(-R), . . . , log(-R)) whose volume is equal to the volume of a simplex with n orthogonal vectors of 
length nlog(.R) — log(fc), i.e., 



(nlog(fl)-log(*)) n 
vol(Vfc) = 



The height of Vk is given by 



Hence 



n 



vol 



nlog(R) — log(fc) 
= 



nvol(Vfc 



h k 

(nlog(22)-log(fc)r 



(n-1)! 



We can finally write 



< y: f (i2) 

S,k>l 

wa k vol{Sk) 

2^ =i vol(A log )^ 1 ; 

= V Wak ^ (14) 

s yi 1 ( n ~ l)!vol(A log )A; s ^ ) 

•(log(iT) - log(fe))"- 1 . (15) 

Let us now develop this sum further by using the Dedekind zeta function. The Dedekind zeta function (k(s) is 
first written as a Dirichelet series defined by the coefficients a*;, 

oo 
k=l 

Using the fact that 

' a k (log(k)) m 



k s 
k=l 
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and the binomial expansion 



\n-l 



(log(iT) - log(*0) r 

( _ )(iog( J R"))"- i - m (i og (fc)r 



m=0 



we can write an upper bound to Sk(s,S) as a function of Ca^ 5 ) an ^ i ts derivatives at s. We denote 

W\/n 



(n - l)!vol(A log ) 



Now 



S, fc=l 



5, k= 

n.-l 



n-1 , 

a* \ ^ I n — 1 



£ 



5, fe=l L "1=0 

n-1 



A- 

n - 1 



m 



(log( J R"))"- 1 -" l (log(fc))' 



(log(iT)r- 1 -™ x 

5, fc=l 



a fe (log(fc)) r 



Next let us upper bound the latter sum in ([Tol by letting the sum go to infinity: 



£ 

S,k=l 



Qfc(log( fc)) 



< ga fc (log(A;)y 
fc=i 



A:* 



= l^ m, C*tol- 

Inserting this to (fT6l ). we finally get 

n_1 / — 1\ 
< K X ( (log^"))"- 1 -" 1 !^^^)!. 

m=0 V m / 

Concentrating on the leading term we further get 

M^^iOog^r- 1 !^)!- 



(16) 



(17) 



Accurate numerical functions enable to easily evaluate this bound. For instance, one can use the Taylor expansion 
of the Dedekind zeta function implemented in Sage in order to calculate the derivatives. Note that the sum only 
contains n terms, which indeed makes the computation feasible. 

VI. Examples of the accuracy of n k 

Let us consider the field extension K/Q with K = Q(y5) of degree n = 2. As K is totally real, we have 
w = 2. The regulator of the field is p K = 0.481211825059603. 

Let us first set R = 10 (see Figure 2), i.e., 1 < k < 100, and denote by b k the exact number of lattice points of 
absolute norm k within the square centered at the origin with the edge length 2R, i.e., 



6 fe = #{xe AnB\N F/Q (x) = k}. 
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Fig. 1. Illustration of the logarithmic lattice with n — 2 and R — 10. 

n_k (circle) vs b_k (triangle) 
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Fig. 2. Illustration of the estimation error, n — 2 and R — 10. 

The values of rife (without the floor operation), and fk(R) = U n fc — kfc|J (the line connecting the previous two) 
are collected in Figure 2. We can see that the error 

f k (R) = [\n k - b k \\ < 2 

for all k. The values are only given for those k for which a k ^ 0, that is, there exists a principal ideal of norm k. 
For all other k we have n k = b k = fk(R) = 0. 
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n_k (circle) vs b_k (triangle) 




nonzero norms 



Fig. 3. Illustration of the estimation error, n = 2 and k > 2000. 

n_k (circle) vs b_k (triangle) 




nonzero norms 



Fig. 4. Illustration of the estimation error, n — 4 and k > 10000. 



When we increase the size of the constellation by considering norms up to k = 2000, i.e., R = 2000, we still 
have 

f k (R) = LK-MJ <3 Vfc, 

see Figure 3. 
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Fig. 5. The exact values take a staricase form. 




Fig. 6. The estimates form a smooth function. 



In Figure 4 we depict the same plot for a lattice Q(a)/Q with dimension n = 4 and minimal polynomial 
f a (x) = x 4 — x 3 — 3x 3 + x + 1 and for norms up to k = 10000. As we can see, the points are well aligned and 
the error remains small, 

f k (R) = [\n k - b k \\ < 3 Vfc. 

In order to see what happens to the size of error in the estimate n k when the dimension grows, let us consider 
a case with n = 8, w = 2. This is already quite a high delay in practice, as we require encoding over eight time 
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Fig. 7. The frequency of estimation errors for k, 1 < k < 65536. Edge length = 2R = 10 and n = 8. 
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Fig. 8. The cumulative frequency of estimation errors for k, 1 < k < 65536. Edge length = 2R — 10 and n = 8. 



instances. 

The construction is based on the maximal totally real subneld of the 32nd cyclotomic field, namely 

^ = Q(C32 + C 32 1 ) 

having a regulator p K = 28.4375954169998. 

While the absolute error gets bigger when dimension grows, it is still negligible considering that out of all the 
cases more than half correspond exactly to the estimate (i.e., nk = bk, meaning no error), and in the rest of the 
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Fig. 9. PEP/c as a function of SNR using 1) the estimate n& and 2) the exact values bk. Extension Q (a) /Q with f a (x) = x 4 — x 3 — 3x :i +x+l, 
n = 4, R = 10, 1 < k < 10000. 

cases (meaning an error occurs) either the error is very small or (a bigger error) occurs very rarely. In Figures [7] 
and [8] we have depicted the frequency and cumulative frequency of errors, respectively. One can see that cumulative 
frequency as high as 90% is achieved already by errors of size < 15. 

Remark 6: It is worth pointing out that the error induced by the estimate is partly due to the fact that we are 
essentially trying to estimate a staircase function by a smooth function, see Figures 5 and 6. 

Finally, we plot the theoretical estimated and exact PEP curves P e {SNR) as a function of SNR. Figure [9] shows 
that there is no penalty by using the estimates in place of the exact values bk when computing ©; the curves 
are perfectly aligned. We have ignored the constant c since this would be the same for both sums, and simply 
plotted P e = £ £i° = °r f|vsP e = i r £fcT F> lettin § 7 take values over an SNR range. 

VII. Conclusions and future work 

In this paper, we have considered number field lattice codes and provided probability bounds for the PEP and for 
the probability of a correct decision of the eavesdropper on a fast fading channel using Dedekind zeta functinos. 
More refined bounds were then derived using geometric analysis on logarithmic lattices. As a byproduct, an estimate 
for the number of constellation points with certain algebraic norm was given, and its accuracy was demonstrated 
through practical examples. 

Future work will consist of improved analysis on the error made in the estimation and of generalizing the results 
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to complex lattices and MIMO channels. One promising approach is offered by division algebras, along the same 
lines as in HI, H3. 
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